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Spatially discrete reaction-diffusion equations with 

discontinuous hysteresis 

Pavel Gurevich* Sergey Tikhomirov^ 

Abstract 

We address the question: Why may reaction-diffusion equations with hysteretic 
nonlinearities become ill-posed and how to amend this? To do so, we discretize the 
spatial variable and obtain a lattice dynamical system with a hysteretic nonlinearity. 
We analyze a new mechanism that leads to appearance of a spatio-temporal pattern 
called rattling-, the solution exhibits a propagation phenomenon different from the 
classical traveling wave, while the hysteretic nonlinearity, loosely speaking, takes a 
different value at every second spatial point, independently of the grid size. Such a 
dynamics indicates how one should redefine hysteresis to make the continuous problem 
well-posed and how the solution will then behave. In the present paper, we develop 
main tools for the analysis of the spatially discrete model and apply them to a prototype 
case. In particular, we prove that the propagation velocity is of order as t —)• oo 

and explicitly find the rate a. 
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1 Introduction 

1.1 Background 

Hysteresis, or, more generally, bistability, refers to a class of nonlinear phenomena which are 
observed in numerous real-world systems. It arises in description of ferromagnetic materials, 
shape-memory alloys, elasto-plastic bodies, as well as many biological, economical, and social 
models, see [8l [2n] - l2^l?rj . The primary goal of the present paper is to analyze a new mecha¬ 
nism (which we call rattling) for pattern formation in spatially discrete systems of reaction- 
diffusion equations (lattice dynamical systems) with hysteresis. The phenomenon occurs in 
any space dimension, including dimension one, and persists even for scalar equations. As 
it is explained below, our results are relevant not only for lattice dynamical systems, but 
also for continuous systems with hysteresis. On the other hand, they link pattern formation 
mechanisms in hysteretic and bistable slow-fast systems. 

Let us begin with the prototype spatially continuous problem 

{ Vr = Vxx + y-iv), a: e (-1,1), r > 0, 

[n(a;,0) = (^(x), a:e(-l,l), 

supplemented with, e.g., Neumann boundary conditions. Here 'H(-) is the simplest hystere¬ 
sis operator, namely, the non-ideal relay or bistable switch, see Fig. II.H a and the (slightly 
modihed) rigorous dehnition in Section [2l 



Figure 1.1: Hysteresis a) with thresholds a < (3. b) with Figure 1.2: Bacteria den- 
thresholds a = — oo and /3 = 0. sity at end of experiment. 

Hysteresis is dehned by two thresholds a < (3 and two values hi,—h 2 € M (in what 
follows, we are interested in the case hi > 0 > —h2). Given a continuous input function 
w{t), its output 'H{w){t) remains constant unless the input achieves the lower threshold a 
or the upper threshold f3. In the former case, the output either switehes to hi if it was equal 
to —^2 “just before” or otherwise remains hi. Analogously, in the latter case, the output 
either switches to —^2 if it was equal to hi “just before” or otherwise remains —h2. Since 
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the function v{x,t) in fll.ip depends not only on r, but also on the spatial variable x, one 
dehnes 'H(n) = 'H{v{x, •))('r) “pointwise”, i.e., for each fixed x. Thus, the hysteresis operator 
H becomes spatially distributed. 

Problem fll.ip is the simplest model of a reaction-diffusion process in which a diffusive sub¬ 
stance with density v{x, r) interacts in a hysteretic way with a non-diffusive substance that 
affects the diffusive one via the reaction term taking values hi or —h 2 . The hrst model of such 
a type was suggested by Hoppensteadt and Jager [TH]. It consisted of two reaction-diffusion 
equations and one ordinary differential equation and described the concentric rings pattern 
that occurs in a colony of bacteria (Salmonella typhimurium) on a Petri plate (Fig. 11 . 21 ) . 

Numerical simulations in [I61IT7] yielded a pattern that was consistent with experiments, 
however the rigorous mathematical description of the model was lacking. To begin with, 
the well-posedness was an open question, due to the discontinuous nature of the hysteresis 
operator. First analytical results were obtained in [21122] (see also [2l[T2l|27| and a recent 
survey [2^), where existence of solutions for multi-valued hysteresis was proved. Formal 
asymptotic expansions of solutions were recently obtained in a special case in [T8| . Questions 
about the uniqueness of solutions and their continuous dependence on the initial data as well 
as a thorough analysis of pattern formation still remained open. 

In [T21[T5] . we formulated the so-called transversality condition for the initial data (p{x) 
in fll.ip that guaranteed existence, uniqueness, and continuous dependence of solutions on 
initial data for scalar equations with hysteresis. In [T3], this condition was generalized to 
systems, and in [101 to the case x G For problem fll.ip . the transversality loosely speaking 
means that if (p{xo) = a or (p(xo) = (3 for some xq G (— 1 , 1 ), then (p'(xo) ^ 0 . Due to [T2IIT4] . 
either the solution exists and is unique for all r G [0, oo), or there is T > 0 such that the 
solution exists and is unique for r G [0,T] and n(x,T) is not transverse. The approach 
of [T 2 HT 3 I was based on treating the problem with transverse initial data as a special free 
boundary problem. The study of regularity of the emerging free boundary was initiated 
in m. For an overview on classical free boundary problems of both elliptic and parabolic 
types, we refer the reader to [SIE112S] and the references therein. 

The key question which we address in this paper is how the solution may behave after 
it becomes nontransverse. To answer this question, we consider the nontransverse initial 
data. First, set /3 = 0 (without loss of generality) and consider an initial function (p(x) = 
—cx^ -f- o(x^) in a neighborhood 13(0) of x = 0. By taking a smaller neighborhood if needed, 
we have (p(x) < 0 for x G B{0) \ {0}. We dehne the hysteresis at the initial moment in this 
neighborhood as follows: = — ^2 and 'H{ip{x)) = hi for x 7 ^ 0. Now we “regularize” 

the parabolic equation in B{0) by discretizing the spatial variable: for any £ > 0, setting 
Vn{T;e) := v{en,T), we replace the continuous model fll.ip in B{0) by the discrete one 


dr 

w(o) 


-f- Hivn), r > 0, n = -N^, ■ ■ ■, iVe, 

-c{En)^ + n ^-N^.... 


( 1 . 2 ) 


where := Vn+i — -|- Vn-i and —)■ cxd as e —)■ 0. Since we are interested in small 

e and in the behavior near the threshold /? = 0 (i.e., in a small neighborhood 13(0)), we 
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consider the next approximation by omitting o{e‘^n^) in the initial data, replacing by oo, 
and formally setting a := — oo. Thus, fll.2p assumes the form 





+ T-L{vn)-, r > 0, n e Z, 


'yn(O) = —c{enY^ n E 7^, 


(1.3) 


the hysteresis operator is represented by Fig. II. li b (see the rigorous dehnition in Section H]). 
A nontrivial dynamics occurs in the case hi > 2c > 0 > —h 2 . To indicate the difficulty, 

note that, due to the initial conhguration of hysteresis, we have -^(0;^) = —h 2 — 2c < 0, 

dr 


but ^-^(0; c) = hi — 2c > 0 for n G Z \ {0}. Thus, for small r > 0, VqIt] e) decreases, while 
dr 

all the other nodes Vn{T;e), n E 7 \ {0}, increase. It is not clear at all, which node achieves 
the threshold (3 = 0 and switches hrst and hence what a further dynamics is. 

In fact, numerical analysis does not reveal any general rule that could describe the behav¬ 
ior of Vn{T] e) for small r. However, it reveals the formation of quite a specihc spatio-temporal 
pattern for large r, see Fig. 11.31 



Figure 1.3: Upper graphs represent spatial prohles of the hysteresis 'H{vn) and lower graphs 
the spatial prohles of the solution a) Nontransverse initial data, b) Spatial prohles at a 
moment r > 0 for h 2 = 0. c) Spatial prohles at a moment r > 0 for h 2 = hi > 0. 

If ^2 = 0, then each node eventually achieves the threshold (3 = 0 and thus 'H^Vn) 
eventually switches from hi to ^2 = 0 for each n E 7. If ^2 > 0, then some nodes achieve 
the threshold and some do not. If we denote by Ni{j) and N 2 {j) the number of nodes in 
the set {no,n±i,... ,v±j} that switch and do not switch, respectively, on the time interval 
[0, cx)), then numerics suggests that 


j^oo Ni{j) 


h2 

hi' 


(1.4) 


Moreover, if h 2 /hi = P 2 /P 1 , where pi and p 2 are co-prime integers, then, for any j large 
enough, the set {uj+i,..., contains exactly pi nodes that switch and p 2 nodes that 

do not switch on the time interval [0, cxd). 

The next numerical observation is as follows. Let = Tn{e) be the switching moment 
of the node n„(r; e) if this node switches on the time interval [0, cxd) and := 00 otherwise. 
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Then, for any fixed h 2 > 0, the r„’s that are hnite satisfy, as n ^ cx). 


Tn 


a{enY + 



if h2 = 0 , 
if ^2 > 0 , 


(1.5) 


where a > 0 depends on hi/c but does not depend on h 2 or £ and O(-) does not depend on £. 

Remark 1.1. In Section [L^ we will show that £ in fll.3l) can be scaled out, see scaling fll.Sp . 
In particular, all the numerical observations concerning the dynamics of have been done 
for £ = 1 and then transferred to an arbitrary e according to the scaling in fll.81) . 

Consider the function 


H{x,T;e) :='H{vn{-;e)){T), x ^ [en — e/2,en + e/2), n E Z, 

which is supposed to approximate the hysteresis l-i{v{x, ■)){t) in fll.ip . Assuming the dy¬ 
namics (11.41) and (ll.hh and taking into account Remark 11.11 we see that H{x,T-,e) has no 
pointwise limit as £ —)■ 0, but converges in a certain weak sense to the function H{x, r) given 
by H{x, r) = 0 for r > and H{x, r) = hi for r < ax^. We emphasize that H{x, r) does 
not depend on h .2 (because a does not). On the other hand, if h 2 > 0, the hysteresis operator 
n{v{ X, •))('^) fll.ip cannot take value 0 by dehnition, which clarihes the essential difficulty 
with the well-posedness of the original problem fll.ip in the nontransverse case. To overcome 
the non-wellposedness, one need to allow the intermediate value 0 for the hysteresis operator. 

Such a re-de£nition of hysteresis is consistent with the behavior of v{x, r) (also observed 
numerically) in the following sense. For a fixed e > 0, the spatial prohle of n„(-;e:)(r) forms 
two humps propagating away from the origin according to fll.5l) . The cavity between the 
humps has a bounded steepness characterized by the relations 

\vk+i{T-,e) - Vk{T-,e)\<he'^, |/c| < n, r > r„, n = 0 , 1 , 2 ,..., ( 1 . 6 ) 

where 6 > 0 does not depend on k, n, and £. As time goes on, the prohle executes down¬ 
wards and upwards motions, always remaining beneath the threshold /3 = 0 and hitting this 
threshold at specihc nodes characterized by fll.4p . We call such a behavior of Vn and H^Vn) 
rattling. Furthermore, numerics indicates that, as £ —)■ 0, the function 

V (x, r; e) := Vnir; e), x E [en — e/2, en + e/2), n E Z, 

approximates a smooth function V (x, r), which satishes V (x, r) = 0 for r > ax^ due to fll.51) 
and fll. 6 l) . In other words, V{x,t) sticks to the threshold line /3 = 0 on the expanding 
interval x E {—i/rja, i/rja). 

We recall paper |3], in which Alt proved the existence of a function V{x, r) that satishes 
the equation 

Vr = I4x + 7(x,'r), 

where 7 (x, r) = ?{(V(x, ■))(t) a.e. on the set A := {(x, r) : V(x, r) ^ a, /d} and 7 (x, r) = 0 
a.e. on the set B := {(x, r) ■. V{x,t) = a or /?} (which potentially may have a nonzero 
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measure). Thus, our heuristic argument provides a qualitative description of the sets A 
and B and justihes the completion of hysteresis by the zero value via the thermodynamical 
limit. To make this argument mathematically rigorous, we should hrst rigorously describe 
the rattling phenomenon in the discrete system fll.Sp . This is the central topic of the present 
paper, in which we concentrate on the case ^2 = 0 and develop general tools for treating 
discrete reaction-diffusion equations with discontinuous hysteresis. The application of these 
tools to the case ^2 > 0 will be a subject of a forthcoming paper. We expect that these tools 
will be applicable whenever h 2 /hi is rational. 

Before we proceed with the description of our tools and of the structure of the paper, 
let us make two more comments. First, the rattling phenomenon also occurs in multidimen¬ 
sional domains. For example. Fig. [T3] illustrates the switching pattern for a two-dimensional 
analog of (11.31) . where we have implemented spatial discretizations on the square and trian¬ 
gular lattices, respectively. Moreover, numerical analysis of the Hoppensteadt-Jager system 
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Figure 1.4: A snapshot for a time moment r > 0 of a two-dimensional spatial prohle of 
hysteresis taking values hi > 4c > 0 and h 2 = hi > 0. The nontrasverse initial data is given 
by = —c(xi + X 2 ). Grey (black) squares or hexagons correspond to the nodes that 
have (not) switched on the time interval [0,r]. a) Discretization on the square lattice, b) 
Discretization on the triangular lattice. 

indicates that the solution remains transverse as long as the central disc in Fig. 11.21 gets 
formed, but the formation of all the rings occurs via rattling. 

Second, the rattling phenomenon is not a pure consequence of a discontinuous nature of 
hysteresis, but rather a consequence of bistability in a system. In particular, it persists in 
bistable slow-fast reaction-diffusion systems. The simplest example is the system 

Vr = V^x + W, 6Wr = f{v,w), (1.7) 

where 5 > 0 is a small parameter and the nullcline of f{v,w) is Z- or S'-shaped. Formally, 
system fll.7p can be treated as another regularization of system fll.ll) . In the case where the 
nullcline of f{v,w) is S'-shaped, one should replace hi and —^2 in the dehnition of hysteresis 
'H(n) by appropriate functions Hi{v) and H 2 {v), see Fig. 11.51 

As 5 —)■ 0, the spatial prohles of v and w in fll.7p behave similarly to V{x,T;e) and 
H{x,T;e), respectively, as £ —>■ 0, see Fig. 11.61 with the exception that the prohle of w 
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remains continuous and forms steep transition layers between mildly sloping steps of width 
tending to 0 as 5 —)• 0. Interestingly, the time-scale separation parameter S in fll.7p yields 
the same effect as the grid-size parameter e in fll.dp . As far as we know, such a rattling 
phenomenon for slow-fast systems has not been explained in the literature, either. 



Figure 1.5: a) The nullcline of the S'-shaped 
nonlinearity f{v,w). b) Hysteresis with non¬ 
constant branches Hiiy) and H 2 {y). 


tends to 0 
as ^ 0 
-T 



Figure 1.6: Upper and lower graphs are 
spatial prohles of the solution w{x,t) and 
v{x,t), respectively, for problem fll.7l) with 
initial data uI ,-=0 = —cx‘^+o{x‘^), w|r=o = hi. 


1.2 Structure of the paper 

Now we come back to the main topic of this paper, namely, discrete system fll.Sp . As it was 
mentioned in Remark fFTl e in (11.31) can be scaled out. Indeed, setting 

t := Un{t) := e~‘^Vn{r] e) (1.8) 


and using the equalities (recall that a = —oo and (3 = 0) 

'hi{Vn){T) = = 'H{Un){t), 

we can rewrite (11.31) as follows: 

{un = ^Un + 'H{Un), t > 0, U G Z, 

I Mn(0) = —cn^, n G Z, 

where ' = d/dt. Problem (11.91) does not involve e, which justihes the fact that Unit) in (11.81) 
does not depend on e. Note that c in (11.91) could be also scaled out replacing Unit), hi and 
—^2 by cUnit), chi and —ch 2 , respectively. We prefer not to do this, in order to keep track 
of what exactly is influenced in our intermediate calculations by the “tangency” constant c. 

From now on, we concentrate on the case h 2 = 0. Due to (11.51) and (11.81) . the asymptotics 
for the switching moment of the node Unit) is expected to be 

tn = an^ + qn, \qn\ < E^/n, ( 1 . 10 ) 


where E > 0 does not depend on n G Z. 
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Our main result fTheorem I3.2p is as follows. Let hi > 2c > 0 and h 2 = 0. Assume that 

finitely many nodes Un(t) (ill) 

switch at moments n = 0,1, • • •, no, satisfying fll.lOp . 

where the constants a = a{hi/c) > 0 and no = no{E) = no{E,hi,c) will be explicitly 
specihed in the main text. Then each node Un{t), n G Z, switches; moreover, the switching 
occurs at a time moment satisfying fll.lOp . 

Since we will provide an explicit formula for the solution Un{t), the fulhllment of hnitely 
many assumptions (ll.lip can be verihed numerically with an arbitrary accuracy for any 
given values of hi and c (see Fig. 11.71) . 

a) 

0 | . ► 

5 10 15 20 25 n 

-0,5- 

_ I 5 . ♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

- 2 - ♦ 

x/TT 


b) 

0 I . —► 

_Qj: 5 10 15 20 25 30 35 n 

-0,2- 

-0,3- 

-0,4- 

-0,5- ♦. 


- 0,6 ♦ 



Figure 1.7: Values of qn /for c = 1/2 and a) hi = 1.5, b) hi = 2.0. 

The paper is organized as follows. In Section O we give dehnitions for the hysteresis oper¬ 
ator and for the solution of problem fll.Op . Next, we formulate the existence and uniqueness 
theorem (Theorem [23]), which includes a representation of the solution M„(t) via the discrete 
Green function Unit). In particular. Theorem 12.51 implies that Unit) = M_„(t), n G Z. 

In Section [3l we formulate our main result (Theorem 13.21) . 

In Section m we formulate three main ingredients for the proof of the main result. 

1. The hrst ingredient is asymptotic formulas for the Green function Unit) and for its time 
derivatives, which were derived in P. 

2. The second ingredient is three equations for hnding the constant a from equation fll.lOp . 
The equivalence of these equations as well as the existence and uniqueness of their root 
a > 0 are proved in Appendix |A| 

3. The third ingredient is the approximation of some singular integrals by Riemann sums 
and corresponding error estimates, which are proved in |15j . 

Sections IHlini and [7| are three key steps in the proof of our main result. The scheme of the 
proof is inductive. Assume we have proved that to,ti,... ,t„_i satisfy fll.lOp for some hxed 
n > no + 1. We fixt\\e hysteresis conhguration, i.e., set Hn := 'H(nn)(tn-i) and consider the 
solution Unit) of the problem 




























(we abuse the notation by using the same letter v as in Section [LT|) . Obviously, n„(t) = Un(t) 
as long as the nodes Vn(t) , Vn+i(t) , Vn+2(t) ,... remain below the threshold (3 = 0. 

The main theorem of Section 1^ iTheorem I5.4p claims that the equation + g„) = 0 

has a root qn satisfying fll.lOp . To prove this, we use an explicit representation of via 
the convolution of with the Green function i/nit) (see fl5.19p i. Then we use asymptotic 
formulas for ?/„(f) (the hrst ingredient from Section H]) and replace the convolution by a 
singular integral (the third ingredient from Section 0]). As a result, we obtain a leading order 
term of order n^, which depends only on a and hi/c, and a remainder of order y/n, which 
also depends on go, <?i, • • •, Qn-i (that are known due to the inductive hypothesis) and on the 
unknown g„. It appears that the coefficient at vanishes due to the choice of a (the second 
ingredient from Section HP . The hard part is to show that the remainder vanishes for some 
Qn satisfying (ll.lOp . This is done by an application of Brouwer’s hxed-point theorem. 

The time moment tn '■= + Qn given by Theorem 15.41 is a candidate for being the 

switching moment of Un{t). To show that it is the switching moment, we have to prove that 
neither of the nodes n„+i(t), n„+ 2 (t),... achieves the value /3 = 0 on the interval (tn-i,^n], 
while Vn{t) achieves it at the moment for the hrst time. This is done in Sections [6] and [3 
In Section El we prove that Vn+iitn) < 0 (Theorem 16.21) . To do so, we estimate the 
gradient Vvnitn) '■= Vn+iitn) — Vnitn) by using the representation of Vvnit) via the gradi¬ 
ent Vynit) of the Green function, applying asymptotic formulas for Vynit) (recall the hrst 
ingredient from Section HP and again replacing the corresponding convolution by an integral 
(recall the third ingredient from Section HP- It appears that the leading order term of order n 
vanishes due to the second ingredient from Section HI Thus, we calculate the next term in the 
asymptotics, which turns out to be —3hi/4 < 0. Hence, n„+i(G) = Vvnitn) < —3hi/8 < 0. 

In Section[7l we hrst show that n„(t) does not achieve the threshold (3 = Ofort E (tn-i, tn) 
fTheorem 17.ip . To do so, we divide the interval (t„_i,t„) into two parts. We prove that the 
function Vn{t) is so small on the hrst interval that it cannot overcome the distance exceeding 
—3hi/8 (the value coming from Theorem 16.21 with n -|- 1 replaced by n). Then we prove 
that Vn(t) is nonnegative on the second interval. Hence, the equation Vn(t) = 0 has a unique 
root, which must be tn- In particular, Vn(t) < 0 for t G (G_i,G). Finally, we show that 
Vvjit) < 0 for all t G (G_i,G] and j > n, which implies that the nodes Vn+iit), Vn+ 2 (t),... 
remain negative for t G (tn-i,tn] fTheorem 17.31) . 

In Section [HI we combine the results from Sections [5l El and [7] and rigorously implement 
the inductive scheme, which completes the proof of the main result, namely. Theorem 13.21 
The crucial role in our main result iTheorem 13.21) is played by the number uq = nQ{E), 
which determines the number of switchings one has to check “by hand” (see fll.lip b The 
number no{E) is determined explicitly by 12 inequalities that must hold for n > no{E). 
Each inequality is referred to as a requirement and is introduced in the text where it is 
used for the hrst time. These 12 requirements contain constants that are also introduced in 
the text where they are used for the hrst time. For reader’s convenience, we have collected 
all those constants in Appendices IB.lHB.3l and the 12 requirements in Appendix IB.41 

The graphs in Fig. 11.81 represent the values of a, E, and no{E) that fulhll assump¬ 
tions fll.llH for c = 1/2 and hi = 1.1,1.2,1.3,..., 2.5. 
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Figure 1.8: Dependence on hi of the values of a, E, and uq^E) that fulhll assumptions fll.lip 
for c = 1/2. a) The values of a are found explicitly for all hi > 1 (see Section belowh 
b), c) The values of E and uq^E) are calculated numerically for hi = 1.1,1.2,1.3,..., 2.5. 


2 Setting of the problem and a proof of its well-posedness 

For a sequence {unjnez of real numbers, we use the notation 

• ^n+1 1 ^n—1 217^ T U^_|_i. 


Let {Mn(t)}nez be real-valued functions dehned for t > 0. We study the problem 

Un = -|- T-i{un), f > 0, n E 7^, (2.1) 

Wn(0) = —cn^, n E Z, (2.2) 

where c > 0 and 'H{w){t), f > 0, is the hysteresis operator dehned for functions w E (^[0, oo) 
such that ta(0) <0 by 


n{w){t) 


hi if w{s) < 0 for all s E [0,f], 

0 if w{s) = 0 for some s E [0, t], 


(2.3) 


where hi > 0 is hxed. In other words, the output of hysteresis is hi unless the input achieves 
the zero threshold; at this moment, the hysteresis switches and since then the output of 
hysteresis remains 0. In the context of problem fl2.1l) - fl2.3p . we will say “a node Unit) 
switches” or “a node n switches” whenever Unit) achieves the value zero for the hrst time. 

Remark 2.1. In the terminology of, e.g., [2U1E7] . the hysteresis operator fl2.3p is a non-ideal 
relay with the thresholds —oo and 0; see Fig. II.li b. 

From now on, we assume throughout that the following condition holds. 

Condition 2.2. hi > 2c > 0. 

We note that the function 'H(u)(f) may have discontinuity (actually, at most one) even 
if u G C°°[0, oo). Therefore, one cannot expect that a solution of problem fl2.1l) - fl2.3|) is 
continuously differentiable on [0, oo). Thus, we dehne a solution as follows. 
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Definition 2.3. We say that a sequence {Mn(^)}nez is a solution of problem fl2.ip - fl2.3p on 
the time interval (0,T), T > 0, if 

1. Un e C[0,T] for all n G Z, 

2. for each t G [0,T], there exists A,a>0 such that sup |un(s)| < Ae"!"'! for all n E'Z, 

se[o,t] 

3. there is a finite sequence 0 = Tq < Ti < ■ ■ ■ < tj = T, J > 1, such that Un G 'O+i) 

for all n G Z and j = 0,..., J — 1, 

4. the equations in fl2.ip hold in (tj, Tj+i) for all n G Z and j = 0,..., J — 1, 

5. Mn(0) = —for all n G Z. 

We say that a sequence {un{t)}nez is a solution of problem fl2.ip - (l2.3p on the time interval 
(0, oo) if it is a solution on (0, T) for all T > 0. 

Remark 2.4. If {Mn(t)}nez is a solution, then, as we have mentioned above, the function 
'H{un){t) has at most one discontinuity point for each fixed n E Z. Hence, the equations 
in fl2.1l) imply that each function Un{t) has at most one discontinuity point on [0, cxo). 

Before we treat existence and uniqueness of a solution, let us introduce one of our main 
tools, namely, the so-called diserete Green funetion 


Vnit) : = 


1 



_ g—2i(l—COS0) 

2(1 — cos 9) 


de, 


t > 0. 


(2.4) 


One can directly check that G O°°[0, cxo) and Unit) solves the problem 


{yo = A|/o -hi, t > 0, 
lyn = A|/„, t > 0, n 7^ 0, 
[ |/„(0) =0, u G Z. 


(2.5) 


Below, we will use the fact that 


yn+i{t) Kijnit), f > 0, n = 0,1,2,..., 


( 2 . 6 ) 


which follows from the formula yn{t) = e ^0„(2f), where In{s) is the modified Bessel function 
of the first kind (see lU Sec. 9.6.19]), and from, e.g., |1]. We will also use the estimate, which 
follows from the series representation of the modified Bessel function PQ Sec. 9.6.10]: 


0 < 2/n(t) = e ^*/n(2t) = e ^ 


4.2m 


- 


m\{m + |n|)! 


< 


2t^|n| 


n E Z. 


m=0 


n ! 


(2.7) 


Below we prove the following existence and uniqueness result. 
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Theorem 2.5. 1. Problem fl2.ip - fl2.3l) has a unique solution {Mn(t)}nez on the time in¬ 

terval (0, oo). 

2. Let tn be the switching moment of the node Unit) if this node switches on the time 
interval [0, cxd) and tn '■= oo otherwise. Then 

9 

CTl 

tn > 7-n e Z. (2.8) 

hi — 2c 

3. Let S{t) be the set of nodes that switch on the time interval [0,t], i.e., 

S{t)-={keZ-H{uk){t) = t)}, (2.9) 

and let |*S'(t)| be the number of elements in S{t). Then S{t) is finite for each t > 0, 
symmetric with respect to the origin, |S'(t)| —)• oo as t ^ oo, and 

Unit) = -cn^ + {hi - 2c)t -hi ^ Vn-kit - 4), t E [0, oo), (2.10) 

k£S{t) 


where we put yn-kit — 4) = 0 for t < t^, 

4- for each n eN, we have t_n = tn and U-nit) = Unit). 

Proof. Step 1. Using fl2.5p . we see that the functions 

^n\t) '.=—cn^ + ihi — 2c)t — hiynif), n E Z, (2-11) 

satisfy the initial condition fl2.2p and the equation in fl2.ip as long as Zn\t) < 0 for all 
n eZ \ {0}. By comparing Zn\t) with the solution 

zfit) = —cn^ + (hi — 2c)t (2-12) 

of problem fl2.1l) - fl2.3p with PLiun) replaced by hi for all n G Z, it is not difficult to see that 

zi^\f) < zfit), t>0,nEZ. (2.13) 

Therefore, the time moment tn at which Znit) vanishes for the hrst time is not less than the 
moment tf = cn^/(hi — 2c) at which zfit) vanishes. 

In particular, Zn\t) < 0 for all n G Z \ {0} and t E [0, c/(hi — 2c)). Let Ti := sup{t > 
0 : Zn\t) < 0 Vn G Z \ {0}}. It follows from fl2.1ip . (14.31) . and (14.7p that ti is hnite. As 
we have seen, xi > c/(hi — 2c) > 0. Furthermore, Zn \t) satisfy the growth condition from 
item[2]of Dehnition l2.3l for t E [0, Xi]. This follows from (12.lip and the fact that \ynit)\ given 
by (12.41) are bounded on any hnite time interval, uniformly with respect to n G Z. Thus, 
Unit) ■= Zn\t) is a solution of problem (I2.ip - (I2.3|) on the time interval (0,xi). 

Let us prove that the solution Unit) is unique on (0, xi). Assume we have another solution 
Unit) on a time interval (0,xi), where xi < xi is such that Unit) < 0 for all t E (0,xi) and 
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n E Z \ {0}. Then the difference Wn(t) := Unif) — Unit) must satisfy the homogeneous 
diffusion equation on the time interval (0,fi) with the zero initial data 

Wn{t) = /\Wn{t), t e (0,fi), neZ, 

'?i'n(0) = 0, n E Z. 

If we looked for solutions that are square summable with respect to n eZ, then the applica¬ 
tion of the discrete Fourier transform would immediately imply that all Wn{t) = 0. However, 
we are interested in solutions that may have exponential growth with respect to n G N (see 
item m in Definition [23]) • We will argue as follows. For each N eN, we consider the functions 

Cn{t) = Cnit) ■= Wn{t) for |n| < N, Cn{t) = (n (t) ■= 0 for |n| > iV 1. (2.14) 


They satisfy the relations 

Cn{t) = AC„(t) + Gn{t), tE (0, Ti), nE Z, 

Cn(0) = 0, HEZ, 

where Gn (t) = 0 for |?7,| < iV — 1 and |?7,| > iV -|- 2, G^^{t) = w±{N+i)it), and = 

—w±N(t). Since no more than finitely many elements in the sequences {Cn(t)}nez and 
{Gn {t)}n£z are nonzero, we can apply the discrete Fourier transform to fl2.15p and obtain 

v+i 

Cnii^)=Yl / Vn-kit-s)G^{s)ds, t E[0,fi], n eZ. (2.16) 

\k\=N''^ 

Now let us £x n > 0 and t E [0, Ti]. By assumption, there exist H, a > 0 such that 

sup |u;fc(s)| < k E Z. (2.17) 

se[o,4] 

Combining fl2.14p . fl2.16p . fl2.7p . (12.17p and choosing N > n, we have 

-^+1 fN aN 

l“n(*)l = IC"(*)I £ * sup (|!/„_i(s)| ■ |Gf(s)|) < c—--- -» 0 as Af -» oo, 

(A'-n)! 

where c = c(n, t) > 0 does not depend on N. Therefore, W-nit) = Wnif) = 0. This proves 
that Unit) is a unique solution of problem fl2.1l) - fl2.3p on the time interval (0, ri). 

Step 2. Set S'(ri) := {0}U{n E Z : Zn\Ti) = 0}, cf. (12.9p . Due to (I2.12p and (12.131) . the 
set S{ti) is finite. Due to (12.111) and the symmetry Unit) = y_n{t), the set S'(ri) is symmetric 
with respect to the origin. Note that to = 0 is the switching moment of the node 0, while 
tn = Ti are the switching moments of the nodes n E S{ti)\ {0}. 

Using (12.51) and assuming Unit) := 0 for t < 0, we see that the functions 


4^Ht) = -cri^ + (hi - 2c)t - hi I i/„(t) ^ Vn-kit - Ti) 

y fceS(Ti)\{o} 

=-CTi^ + {hi - 2c)t - hi ^ yn-k{t-tk), tE[0,oo), 

k£S{Ti) 


(2.18) 
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satisfy the equations in fl2.ip as long as Zn\t) < 0 for all n G Z \ S'(ri), i.e., as long as 
S{t) = S{ti). Obviously, Zn\t) also satisfy the initial condition fl2.2p . 

As in Step 1, we see that the time moment at which Zn{t), n G Z \ S'(ri), vanishes 
for the hrst time is not less than cn^/{hi — 2c). Hence, there is a positive time interval (of 
length bigger than Ti) on which Zn \t) < 0 for all n G Z \ 

Let T 2 := sup{t > 0 : Zn\t) < 0 Vn G Z \ It follows from fl2.11l) . fl4.3p . and fl4.7l) 

that T 2 is hnite. We have proved that T 2 > ti. Furthermore, Zn\t) satisfy the growth 
condition from item [2] of Dehnition 12.31 for t G [0,r2]. Thus, Un{t) := Zn\t) is a solution of 
problem fl2.ip - fl2.3l) on the time interval (0,r2). Note that Un{t) = Zn\t) for t G [0,ri]. The 
uniqueness of Un(t) on the interval {ti,T 2 ) can be proved similarly to Step 1. 

Continuing these steps, we obtain the desired inhnite sequence {Tj}j>o from Dehnition l2.31 
On each step, we compare Unit) with z'^it) given by (I2.12p and conclude that the switching 
moments satisfy > cu?/{hi — 2c). Hence, Tj —)■ oo as j —>■ oo. □ 

3 Main result 

We recall that Condition 12.21 is assumed to hold throughout. Below in the text we dehne 
a > 0 (see Lemma . Eq > 0 (see fl5.12p i and an increasing function no : {Eq,oo) —?• N 
(see Requirements [THT2] in Section IB.4P . 

Definition 3.1. We say that a number E > Eq is admissible if the following holds: 

1. each node u^, k = 0, ±1 ..., ±no = ±no(i?) switches at a moment tj. satisfying 

tk = ak^ + qk, \qk\<Ey/n^, (3.1) 

while neither of the nodes M±(no+i))'W±(no+ 2 ), ■ ■ • switches on the time interval 

2. at the switching moment tn^, we have 


^no+l(^no) V^no(^no) — 


3hi 


(3.2) 


The main result of this paper is as follows. If hnitely many nodes fc = 0,..., nQ{E) switch 
at time moments satisfying fl3.ip . then all the nodes n G Z will switch and their switching 
moments will be of order an^. On Fig. 11.81 b. one can see the values of admissible E, which 
we found numerically for c = 1/2 and hi = 1.1,1.2,1.3,..., 2.5. Figures fL8l a and 11.81 c 
depict corresponding values of a and no(£'), respectively. 

The rigourous formulation of our main result is as follows. 

Theorem 3.2. Assume that E > Eq is an admissible number, and let uq = nQ{E). Then 
for all n > Hq + 1: 
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(3.3) 


1. Each of the nodes Un switches at a moment tn satisfying 

tn = an^ + qn, \qn\ < 
tk ^ tno ^ ^no+1 ; k 0, 1, , 77-0 1, 


2. There exists > 0 depending on hi,c,E, but not on n, such that 


^Unitfi) E 


3/ii 


< A^n VUn{tn) < 


3/ii 


4 Auxiliary Statements 

In this section, we formulate several auxiliary statements. Each of them is a key ingredient 
in the proof of our main result, i.e.. Theorem 13.21 

In Section ITT] fProposition [4. ip . we establish asymptotic formulas for the discrete Green 
function yn{t) given by (12.dh . It is essential that the leading order terms in the asymptotics 
depend only on n/\/f, while the remainders are estimated uniformly with respect to n. 

In Section 14.21 we consider three expressions containing integrals fl4.12p of leading order 
terms in the asymptotics of ynif), Vyn{t), and yn{t), respectively. These three expressions 
will enter the leading order terms in asymptotic formulas for Unifn), Vunitn), and iinitn)- In 
Proposition 14.21 we show that these terms vanish for the same value of a, thus determining 
the “propagation rate” an^ in the switching moment asymptotics for tn in fl3.ip and fl3.3p . 

In Sectionwe elaborate on properties of integrals fl4.12p from Section |T]2l In the proof 
of our main result, these integrals will play the role of approximation of some Riemann sums. 
Note that the corresponding integrands are not smooth functions, but have singularities of 
order (1 — or (1 — x)~^k^ at a; = 1. In Propositions 14.31 and 14.41 we provide error 
estimates for approximation of such integrals by their Riemann sums. 


4.1 Properties of the discrete Green function yn{t) 

Consider the functions h, f,g,f: M+ —)■ M given by 


4, f{x):=2xjy^h{y)dy, g{x) := f{x), f{x):= 


h'”{x) 

6x 


(4.1) 


Note that these functions belong to G°°[0, oo) and decay to zero as x —)■ cx), together with 
all their derivatives, faster than any exponential. Moreover, 

h{x) = g’{x) = f{x), 2h{x)+xg{x) - f{x) = 0, ^(0) = -^. (4.2) 

Consider the functions rQ,ri,ri,r 2 ,wo,wi : (M U {0}) x M+ —)■ M satisfying the following 
relations for n = 0,1, 2,... and f > 0: 



+ ro{n,t), 



+ ri{n,t), 


(4.3) 
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yn{t) = +ri(n,t), Vnit) = +r 2 (n,t), (4.4) 

V?/n(t) = +wi{n,t), Vynit) = g\^j + J +wo{n,t). (4.5) 

We fix throughout the paper 

To > 0. (4.6) 

The following estimates are proved in Un- 

Proposition 4.1. There exist constants Aq, Ai, A 2 , Ai, Bo, Bi, A2, B2 > 0 {depending on tq) 
such that, for all t > tq, n = 0 , 1 , 2 ,, and i = 0 , 1 , 2 , the following inequalities hold: 


1 

\ri{n,t)\ < 

\ri{n,t)\ < 

(4.7) 

\wQ{n,t)\ < Boj, 

\wi{n,t)\ < 

(4.8) 

|ij„(()| < 

Vf/n(t) < i?2^- 

(4.9) 


4.2 Equivalence of some equations 

Consider the functions 


F{a,x) := \/a{l - x"^) f 
1 


H{a, x) := 


y/a{l - a;2) 



where a > 0, a; G (—1,1) and /, g, h are given by (14.Ih . Set 


Ipia) := / F{a,x)dx, Ida) := / G{a,x)dx, Inia) := / H{a,x)dx 


(4.10) 

(4.11) 


(4.12) 


J-1 J-1 J-1 

The following proposition is proved in Appendix lAl 
Proposition 4.2. Each of the three equations 

—c + (hi — 2c)a — hilp{a) = 0, (4-13) 

-2c-hilda) = 0, (4.14) 

(hi — 2c) — hilnia) = 0 (4-15) 


has a unique root on the interval (0, cxd). Moreover, all these equations have the same root. 

In what follows, we £x a given by Proposition 14.21 and write F{x), G{x), H{x), Hi{x), 
If. Ig, Ih, omitting the dependence on a. 
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4.3 Error estimates for Riemann sums 

Let iV G N, and let 2:0 = or = 0. The following propositions are proved in na. 
Proposition 4.3. Assume that a function Fi{x) can he represented as 

Fi{x) = ci(l - xY^'^ + 02(1 - xY^‘^ + Fi{x), 

where Ci,C 2 G M and Fi G C‘^[zo, 1]. Denote the error estimate of the Riemann sum of the 
integral f^^Fi(x)dx by 

“^^0 \ k=ZQn+l ^ / 

1. There exists Li = Li{Fi,N) > 0 such that 

\Rn\ < -^ 1 ^ 7 ^, n> N. 

2. If, additionally, Fi{x) = 03(1 — + 04(1 — + Fi{x), where 03,04 G M and 

F G CYzo, 1], then there exists Li = Li{Fi,N) > 0 such that 

\{n + lYRn+i-n^Rn\ <Li^j- n>N. 

I I 


Proposition 4.4. Assume that a function F 2 {x) can he represented as 

F2{x) = oi(1 - + 02(1 - xY^‘^ + F2{x), (4.16) 

where oi > 0, 02 G M, and F 2 G C^[zq,1]. Then there exists L 2 = L 2 {F 2 ,N) > 0, L 2 = 
L 2 {F 2 , N) > Cl and I 2 = 12 {F 2 ,N) such that 


l 2 -< 

n 


^ 21/2 


F 2 {x)dx 


' ZQ 


n \n 

k=zon 


< L 2 - 


n 


1 

172 ’ 


n > N. 


(4.17) 


In particular. 


n—1 

E iMS 


k=zon 


is hounded. 


Proposition 4.5. Assume that F 3 G C°[—1,1) and |(1 — xY^^F^ (x)| is bounded on [— 1 , 1 ). 
Then there exists L 3 = L^F^, N) > 0 such that 





n > N. 
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5 Asymptotics for Un{tr 

5.1 Preliminaries 

Set {h is defined in (Id.lh ) 

ha{x) := , Da ;= inf ha{x) > 

Vva/ \Xy/aJ a:e(0,l] 


p := snp 


Kix)x 

x£(0,l] ha^X^ 

Lemma 5.1. 1. 0 < p < 2e“^ < 1. 


$ (a;) ;= a; > 0. 


(5.1) 

(5.2) 


2. The function $p(a;) is nonincreasing on (0,1). 

Proof. 1. Let /i := l/(4a). Then 
2^ha{x) = (l + 

2^/Tth'a{x)x = e~^^^2p, (^—x^ ^ e~^^^2p{x~‘^ — x‘^)e~^^^ '^-x'^) < 2e“^e 

Hence, h'a{x)x < 2e~^ha{x) and p < 2e“L The inequality p > 0 is obvious. 

2. Relations fl5.2p imply the following for x G (0,1]: 


-ip-px* 


= 


h'a{x)x^ — ha{x)pX- 

X2p 


.p-1 


j.p+1 


{h'^{x)x -pha{x)) < 0. 


Hence, $p is nonincreasing on (0,1). 
For any iV G N, set 


1-P 


Dpi'.— 2 2 —1, iV I (1 + — 


1+p 
1\ 2 


1 . 


□ 


(5.3) 


Lemma 5.2. For any N ^ N we have for A; > 1, n > the following inequalities hold 

1 


1 + 


1-P 

1\ — 


k 


Dpi^, 


< Dp2, 

Dpi < - 


1+P 

1\ — 


>l + D 


p2 ) 

n 


Proof. Fix a G (0,1). Note that 

(1 + xY~°' < 1 + x{l — a), X E (0,1]. 

Consider the function Pa{x) := ((1 + a;)“ — l)/x. Inequality (15.611 implies that 


(5.4) 

(5.5) 

(5.6) 


PL(x) = 


a{l + a;)“ ^x — ((1 + x)" — 1) 


< 0. 


x^ 
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Hence, Pa{x) is nonincreasing and, for any ki > k 2 , the following holds: 

Pa{l/h) > Pa{l/k2). (5.7) 

Inequalities fl5.4p are straightforward consequences of fl5.7p for a = (1 — p)/2 and a = 
(1 +p)/2, respectively. Obviously, Dp 2 > 2 ^+^)/^ — 1 > Dpi and Dpi < \/2 — 1 < 1/2. □ 

Set {Dpi, Dp 2 are given by fl5.3p i 

X := Dp2 — Dpi — 2{Dpi + Dp2)Dpi — . (5.8) 

Lemma 15.11 and relations (15.311 imply that x > 0 for large enough N. In what follows, we £x 

77 G N such that x > 0. (5.9) 


Assume that no = no{E) satisfies the following. (We remind that the complete list of 
requirements determining no{E) is given in Section fB.dl ) 

Requirement 1. no > N. 

Set 

Cn ■= —cn^ + (hi — 2c)an^ — hi yn-k{(^{n^ — h^)). (5.10) 

\k\<n—l 

In other words, the values Cn are obtained by formally substituting t = an^, tk = ak'^, and 
S{t) = { —(n — 1),..., n — 1} in (I2.inp . In Section [AS] below, we will prove the following. 

Proposition 5.3. There exist K,K' > 0 such that for n > N the following inequalities hold: 


1 


\Cn+l-Cn\<K — , \Cn\<K'V^ 


n 


(5.11) 


Fix 77 and K' from Proposition 15.3! and set 

K + xK' 


Eo := 


° (hi — 2c)x 
Note that Eo> t) due to (15.9p . For each E > Eq, set 


(5.12) 


off™ := a 


2^\/^ I 2^x/n 

2n-l’ •““+2n-l- 


(5.13) 


We assume that Hq = no{E) satishes the following requirement. 

Requirement 2. a™™ > ro/(2n — 1) for all n > Hq, where Tq was fixed in (14.61) . 
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Note that Requirement [2] implies that 

a(n^ — k'^) — 2E\/n > tq, n > hq, \k\ < n — 1. (5-14) 


Set 


2Ri 1 2AlEa^l’^ 1 

aDaU ^ i:)a(a“™)3/2 


(5.15) 


Obviously, 

a™™ —)■ a, —)■ 0 as n —)■ oo, 

a:)7i>ar, <5n+i<5n, n>no. 

We assume that no = no(i?) satisfies the following requirement. 

Requirement 3. + 6 n+i < 2Dpi for all n > Uq. 

Below we will use the following constants Sa, Ta, and R. For a > 0, let Sa be the smallest 
number satisfying the inequalities 


\k\<n—l 


(n — A;)" 

(n2 — A;2^a+l/2 


< Sa—^, n > N. 
Vn 


(5.16) 


For a > 1, let Tq, be the smallest number satisfying the inequalities 

< Tq, —, n > N. 
n“ 


(n2 — k'^) 


Let R be the smallest number satisfying the inequalities 

1 


E 

|fc|<n—1 


(^2 _ p ) l /2 


< i? (= vr), n > N. 


(5.17) 


(5.18) 


5.2 Candidates for switching moments 


5.2.1 Formulation of a theorem on existence of the candidates 

In this section, we will prove the following result. 

Theorem 5.4. Let the assumptions of Theorem \3.2\ hold. Then there exists a sequence t^, 
k E h, such that R, t±i ,..., t±no ore given by fl3.ip . 

tn = an‘^ + qn, \qn\ < E^/n, t_n = R forn>no + l, 

I'k ^ ^ ^ . . . , k 0, 1, ... , TIq 1, 

and the functions 

Vn{t) :=-cn‘^ + {hi - 2c)t - hi ^ yn-k{t-tk), n > Uq + 1, (5.19) 

|/c|<n.—1 


satisfy 


Vn{tn) =0, n > no + 1. 


(5.20) 
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Remark 5.5. The sequence in Theorem 15.41 is a sequence of candidates for switching 
moments in the following sense. Assume that, for some n > no + 1, we know the following 
(this is what we will in particular prove in Sections Eland [7] below): 

1. the nodes uq, ..., Un-i switch at time moments to,, tn-i, respectively, 

2. the nodes Un, Wn+i, Mn+ 2 ,... do not switch on the time interval [0, tn). 

Then Vn(t) coincides with the solution Un(t) of problem fl2.ip - fl2.3l) on the time interval [0, tn) 
and equality (15.201) implies that tn is the switching moment of Unit). 

5.2.2 Proof of Theorem 15.41 

First, we substitute fl5.19p into (15.201) . replace 4 and tn by ak'^ + qk and an^ + Qn, respectively, 
and expand Un-k into the Taylor series around a{n? — k'^). This yields 

0 = - cn"^ + {hi- 2 c){an‘^ + Qn) - hi ^ yn-k{a{n^ - k"^)) 

\k\<n—\ 

-hi ^ yn-k{a{n^ - k‘^)){qn - Qk) 

\k\<n—l 

-hi ]^yn-k{a{n^ - k^) + U,k){qn - qk)'^, n>no + l, (5.21) 

|A:|<n—1 

where < \qn — qk\- We introduce the notation 

0^n,k ■ yn—k{l^{ll k )), (3n,k{qn) ■ '^yn—k{(l{u ) T ^n,k){qn ?fc); (5.22) 

where we omit an explicit indication of the dependence of I3n,k on with \k\ < n — 1 . 
Further, set for fc = l,2, ...,n — 1 

Jn,k{qn) ■ dn,k ~i~ (^n,—k T (iin,k{qn) T Pn,—k{qn)y Jn,o{qn) ■ ®n,0 T f^n,o{qn) j (5.23) 

n—1 

Jn{qn) ■ ^ ^ Jn,ki Dn{qn) ■ hi 2c hiJn{qn). 

k =0 

Using this notation and recalling the dehnition of the constants Cn in (I5.10p . we rewrite 
(I5.2ip as follows (it will also be convenient to replace n by n + 1): 

n 

Cn +1 + Dn+l{qn+l)qn+l + hi E Jn+\,k {qn+i)qk = 0, n> no. (5.25) 

k =0 

Thus, it remains to find a sequence qk, A: G Z, such that \qk\ < Ey/no for fc = 0, ±1,..., ±no, 
|g„+i| < E\/n + 1, g-(n+i) = qn+i for n = no, no + 1,..., and the equalities (I5.25p hold. 
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First, we note that , q±no are already prescribed by the assumption of the theorem. 

Moreover, fl5.25p holds with n + 1 replaced by Uq: 

riQ-l 

Cno + D noiQno )<?no + hi E ■Jno,k {qno)qk = 0. (5.26) 

k=0 

Indeed, Requirement [2] implies that > t^, k = 0,... ,no — 1. Therefore, for all t G 
[tno-i^tno) holds S{t) = {-(no - 1),...,no - 1}, in (I2.10p and 


^no(^'^io) ^^0 (^1 2c)t72o ^ ^ Vno—kitno ^/c)* (5.27) 

|fc|<n.o-l 

Hence, fl5.26p is obtained in the same way as fl5.25p from fl5.19p and fl5.20p . 

Now we proceed by induction. Fix n > uq. Suppose, we have constructed the desired 
sequence qo, ■ ■ ■ ,qn- Let us hnd satisfying |gn+i| < E^/n + 1 and equation (I5.25p . We 
rewrite equation fl5.25l) in the form 


qn+l F(?n+l)) 


F{qn+i) := 


a 


n+l 


Dn +1 {qn+l) 


n 


+J2 


Jn+l,k {qn+l) 

Dji+l {qn+l) 


■ qk- 


(5.28) 


To prove Theorem 15.41 it now suffices to show that if qk G [—E^/n, Ey/n\ for k = 
0, ±1,..., ±n, then F has a hxed point on the interval [—E\/n + 1, E\/n + 1]. 

To do so, we need to show that F maps the interval [—Ey/n + 1, —Ey/n + 1] into itself. 
Let us indicate the main difficulty on this way. We will see in Sections 15.31 and 15.41 that 
Cn ~ y/n, Dn ~ l/y/n, and Jn+i,k{qn+i) ~ l/Vn^ - P, provided that \qn+i\ < Ey/n + 1. 
Therefore, the straightforward attempt to estimate |F(g„+i)| would yield 


|F(Wi)l < 


a 


n+l 


Dn+l{qn+l) 


+ Ey/n ■ hi 

k=0 


Jn+l ,k{qn+l) 
Dn +1 {qn+l) 


(5.29) 


and we would obtain nothing better than |F(g„+i)| < const • n. 

To overcome this difficulty, we will use the following trick. Note that, by the induction 
hypothesis, fl5.25p holds with n + l replaced by n. Therefore, we can multiply fl5.25p by 
1 + x/n with an appropriate x > 0 and subtract fl5.25p with n + l replaced by n. As a 
result, we will obtain the equation 

= F(«„+,), F(,„+i) = - - A, ■ ft, (5.30) 

Dn+l{qn+l) ^^0 Dn+l{qn+l) 


which is equivalent to fl5.28p . The advantage of this new representation will be that we 
will obtain Cn ~ l/\/^ and Dn ~ 1/y/n. Hence, the hrst term in the formula for F 
can be estimated by a constant ai > 0. Furthermore, we will show that the expression 

Jn+l,k{qn+l) 


hi E 

k=0 


Dn +1 {qn+l) 


is estimated by 1 — a 2 ly/n with a 2 > 0. Therefore, fl5.30p will yield 


F(gn+i)| < «! + Ey/n - Ea 2 < Ey/n + 1, 


(5.31) 
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ii E > Eq = ai/a2- In particular, it will turn out that the appropriate x is given by fIS.Sp 
and Eq by fl5.12p . Interestingly, x = 0 would not be sufficient for this scheme as it would 
then follow that 02 = 0. 

To make the above argument rigorous, we need the following proposition, in which we 
do not explicitly indicate the dependence of the functions on qn+i- 

Proposition 5.6. Let the assumptions of Theorem 13.21 hold. Then for any n > uq and 
go, g±i, • • •, Q±{n+i) £ [—E\/n + 1, E\/n + 1], the following holds with x given by (I5.8p : 

L Jn,k> 0, k = 0,1 ,... ,n - 1, 

2 . Jn,k - {1 + ^ Jn+i,k >0, k = 0 ,l,...,n-l, 

3. Dn — hi ^ Jn+l,n ^ 0. 

Now, assuming that Proposition 15.61 is true, we complete the proof of Theorem 15.41 After 
that, in Section l5^ we prove Proposition 15.61 

Substituting n instead of n + 1 into fl5.25p and using fl5.26p for n = no or the induction 
hypothesis for n > uq (and omitting the dependence of Jn,k on g„), we obtain the equation 

n—1 

Cji T EjiQn T h\ ^ ^ Jn,kQk 0. 

A:=0 

Multiplying f|5.25p by 1 + x/n and subtracting the latter expression we have 

0 = {Cn+l — Cn) H- Cn+l + (l + ^ Dn+lQu+l 

n \ nJ 

n—1 

T T k\Jn^l^n L)^ Qn T h\ ^ ^ T Jn+l,k Jn,^ Qk- 

k=0 

Equivalently (cf. (I5.3np ). g„+i = F(g„+i), where F(g„+i) satishes 

^ Dn+lF{qn+l) = — ({Cn+l ~ C'n) H- Cn+lj + (On — ^ hiJn+l,n^ qn 

n—1 

T hi ^ ^ (jn,k Jn+1,}^ qk- 

k=0 

According to Proposition 15.61 all the coefficients at qk, qn, qn+i are positive. The inductive 
hypothesis \qk\ < E^/n, k = 0,... ,n. Proposition 15.31 and the inequality E > Eq imply that 


X 


n 


X 


1 H-) -Dn+l|F(gn+l)| < {Cn+l ~ Cn) H- Cn+l 


n 


n—1 


k=0 

n—1 


Dn Ey/n ^ ^ yJn,k J Jn-\-l,kj Ey/fi 

K + kK' 


n—1 


< 


—= -h E\fn ( (-Cn + hi^ ^ ^n,fc i~hi^H-^ ( Jn+l,n + ^ ^ A 

\\ k=o J n/ 


^n-\-l,k 
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Combining this with fl5.24p yields (cf. fl5.3ip i 


X 


1 + - 


n 


Dn+i\F{q^+i)\ < {K + xK')^ 

Vn 

= {K + xK')^ 
\/n 


+ E^/n (^{hi — 2c) — ^ {{hi — 2c) — Dn+i)'j 

+ Ey/nDn+i H-) — E^/n{hi — 2c) — . 

\ nJ n 


The latter estimate and the inequality E > Eq, where Eo is given by fl5.12p . imply 


|F(q'„+i)I < E^/nDn^i(l + ^ . 

Hence, |F(g„+i)| < Ey/n < Ey/n + 1. Therefore, F maps [—E\/n + 1, Ey/n + 1] into itself. 
Furthermore, the function F is continuous, because the functions 


f3m,k {Qt! 


y^-k{a(m?-k'^)+qm-qk)-y-m-k{a(m?-k'^)) 


Qm—Qk 




0 , 


Qm 7^ Qki 
Qm 


are continuous with respect to G M for fc = 0,..., m — 1 and Dn+i{qn+i) > 0 by Propo¬ 
sition 15.61 Hence, by Brouwer’s hxed-point theorem, the map F has a hxed point on the 
interval [—Ey/n + 1, Ey/n + 1]. The latter implies Theorem 15.41 

It remains to prove Propositions 15.31 and 15.61 which we do in the next two sections. 


5.3 Proof of Proposition 15.3 

Taking into account equalities (14. 3 p and (I4.13p we write (I5.10p as follows: 

— C^(l) C^( 2 ) ^( 3 ) 

~ '-'n + '-'n + 5 


where 



:= hi j n^Ip - 

■ Vain- 


V 

|/c|<n—1 

C7' 

■.= -hi 

|/i;|<n—1 

1 

y/a(n2 — /c2) 

cm 

1 

M 

ri{n — /c, a{n 


\k\<n—l 



n — k 


f\ 


yy/a{n‘^ — k"^)^ 

n — k 


Proof of the second inequality in (15.lip . Below we separately estimate Cn'^ and Cn'^ -|- Cn^. 
Step 1. Proposition 14.31 (item 1) applied to the function E{x) (see fl4.10p h implies that 


IC^I < hiLi{E,N)^=K[^, K[ := hiLi{E,N). 


(5.32) 
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Step 2. Due to fl4.3p . fl4.3p . and fl4.7p . 


|Cf+cf| = 

where, by flS.lSp . 


hi rniji — k,a{n^ — k'^)) 

\k\<n—l 




< ^ 1 , 


li :— hiAo—^R. 


Hence inequality flS.lip is satisfied for 


K' := K[ + h 


y/N' 


(5.33) 


(5.34) 

□ 

Proof of the first inequality in fl5.1ip . Below we separately estimate — \ C^+]^ — Cn\ 

and C^+i — Cn'^ in Steps 1,2, and 3 respectively. 

Step 1. Applying Proposition 14.31 litem 2) to the function F{x), we conclude that 

41) 


Id'i - C“| < A'l := ftiA,(A, iV). 


Step 2. Set 


F{x) := 


--f 


1 — a; 


y/a{l — x‘^) I A/a(l — x‘^) 


, Ip := F{x)dx. 


(5.35) 


(5.36) 


'-1 


Note that limF(a;)(l — x)^/^ < 0. Therefore, by Proposition 14.41 applied to the function 

x^l 

—F{x), for some constants Lj, L*p Ij >0, we have 


L)^ - Ij- < -Ip - -ic'f < L ~. 
t ^/n ^ n hi 1 


Hence, 


ld4 - C<^'l < A24=, I<2 := A, [L,-L) + l, 


n 


|Cf I < Alii ^ 


/ / /^, 

Step 3. Inequalities fl4.7p and fl5.17p imply that 

1 /y 1 ^ r ^ 

a[n^ - P))3/2 - %3/2^3/2^3/2- 

\<n " 


(5.37) 

(5.38) 


Hence, 


|CSi-CW|<A-s^, A3:=2Aiii^r, 


Summarising fl5.35p . fl5.38p . and fl5.39p yields 


a3/2 3/2 
1 


\Cn+l - Cn\ < K — , K-=Ki + K2 + Kh- 
/n N 


(5.39) 

(5.40) 

□ 
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5.4 Proof of Proposition 15.6 

5.4.1 Proof of Proposition [5T6l Preliminaries 

For = — 1, consider the following representation of Jn,k (see fl5.23|) and fl4.4p h 

Jn,k = + {Wn,k + Uln-k) + {f^n,k + f^n-k). 


'^n,k ■ '^n,k^'^n,—ki '^n,±k 

For fc = 0, we have 




ri(n=FA:, {a{n^-k‘^)). 


T _ Tinain , , o rmain__ 

^n,0 ^n,0 ' ^n,0 “T Pnfli '-'n,0 ' ^rifl ■ 


\/i 


-.h 


an^ 


, Wnfi ■■= ri{n,an‘^). 


The general idea is to prove each assertion of Proposition 15.61 for hrst, and then 

consider Jn,k as a small pertnrbation of We formulate this fact as a lemma. 

Lemma 5.7. Let 6n be given by fl5.15p . Then 


IJ, 


n,k Jn^k 
Tmain 


Tmain I 

1 . \ _ n 


< 


A; e {0,l,...,n-1}. 


n + k n — k' 

Proof. Fix /c G {1,..., n — 1} (the case /c = 0 is similar). Due to fl4.7p and fl5.14p . 

1 


(a(n2-P))3/2' 


Dehnition fl5.22p of I3n,k, (I5.14p . fl4.9p and fl5.13p imply that 

\/3n,k + Pn-k\ < , max \yn-k{a{n‘^ - k‘^) + • |g„ - qk\ 

|^n,fc \ ^\Qn~Qk I 


< 


2A2E^/n 


< 


2A*2E^/n 


{a{rp — fc2) — 2E^/n)^P (a™°(n2 — k?)YP 


Note that, by fl5.ip . 


Tmain 




(a(n2-fc2))V2- 


Relations (I5.42p - (l5.44p imply 


|’?^n,fe T ICri,—fc| ^ {a{n‘^—k^)p/^ 2Ai 1 1 


Tmain 

•^n.k 


Da 

(a(n 2 —A: 2 )) 1/2 


aDa n + k n — k' 


2^2 Epn , 

|/ 3 n,A: A l^n,—k\ ^ 2A2EQ, \/n 1 


Tmain 

’^n,k 


Da 


(a(n2—fc2pl/2 


Da(a““)3/2n + fcn- fc' 


(5.41) 


(5.42) 


(5.43) 


(5.44) 


Hence, inequality fl5.4ip holds. 


□ 
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5.4.2 Proof of Proposition 15.61 Part 1 

Since > 0, it suffices to show that the right-hand side in (15.411) is less than or equal 

to 1. The latter is true because < 1 due to fl5.5p and Requirement |3l 


5.4.3 Proof of Proposition 15.61 Part 2 

Fix fc G {1,..., n — 1} (the case fc = 0 can be treated similarly). Note that 


rmain _ 


y/a(n2 - F) \y n + k 


n — k 




n — k 


■Op ( 


n + k j y/a{n — k) {n + k) (i+p)/2 ’ 


where 4’p(a:), x G (0,1], is given by fl5.2p . Taking into account that 4’p(a;) is nonincreasing 
(see Lemma [5Tl item 2) and using fl5.4p . we have 


^ (n + 1 - k)(^-py^{n + 1 + - (i+ ^ 


Tmain 

^n+l,/c 


(n — fc)d-P)/2(72 -|- /^)(l+p)/2 


> I 1 + L) 


pi; 


l + D. 


"n — k J \ ' + k 

Combining fl5.45p with fl5.4ip yields 


>i + n 


n — k 
1 


(l-p)/2 


1 + 


n + k 


(l+P)/2 


pi 


n — k 


+ n 


P 2 


n + k 


(5.45) 


J, 


n,k 


Tmain 

^n,k 


1 + 


J, 


> 


Tmain 

^n,k 


— Tmain 


■I 




1 7 , 1 . — /mam 

n+l,k Jn+l,k 1 + 

I /main 

"'n+l,A: 

5. 


> 


Tmain i _ On n 

^n,k n—k n-\-k 


J: 


1 - 


n 


n — k n + k 


n-\-l,k 


1 + 


1 + 


n-\-l—k n+l+fc 
-1 


^n+l ^ 

n — k n + k 


> ( 1 -|- Dpi ^ + Dp2 


n + k 


1 (*^n T *^n-|-l) 


n 


n + k n — k 


It is easy to show that Requirenment [3] implies that the minimum of the last expression is 
achieved for A; = 0. Hence, using Requirenment [3] again, we obtain (x is given by fl5.8p ) 


J, 


J, 


— > + {Dpi + Zlp2)—^ ~ (^n + <^n+l)—^ 

^ 1 + (^Dpi + Dp2 — {Sn + (5„+l) — {Dpi + Dp2){6n + 5„+l)—^ > 1 + 


X 

n 


5.4.4 Proof of Proposition 15.61 Part 3 

Using (I5.24P and the notation in the beginning of Section 15.4. H we write 

Dn hiJji-ll^n T {hi 2c) hiJfi hiJn+l^n T ^ ,1 (S 2 ^n “t” T U 4 jj), 

^ ^ (5.46) 
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where 


n —1 


:= ih - 2c) JnP. S2,n := hi7n+l,n (l + ^ > 

k=0 

n—1 

Ss,„ := h, Y, (J„,j - J“r) , Ei,„ := fc, + 7„+i,_„) (l + ^) • 

k=0 

In steps 1-4 below, we will estimate ..., We will see that and S 2 ,n are 
“large” with respect to and which motivates the splitting in fl5.46p . Set 

Hi{x) := H{x) + H{—x), x G [0,1), 

where H{x) is given by fl4.1ip . Then for fc = 1,..., n — 1 we have 

JT = -«■ (-). 47” = 

’ n \nj In 

Step 1. Proposition 14.41 applied to the fnnction H{x) (see 04.111) 1 and eqnation 04.12p imply 

> hi (ch—j= — Ih^ , 

V Vn nj 

Ch := Ll{H, N), Ih := kiH, N). (5.47) 

Step 2. Since the function h{x) is decreasing, 

\/2a vn V n/ 

Step 3. Using 05.4ip and Proposition 14.51 applied to the function - -H{x), we have 


1 — x^ 


n —1 

S3,n| < hi6n V < hi5nCH2^, 

^Wn 

k=0 ^ 


where 


Ch2 L?. 


^H2 ■— -^3 1 


I — x^ 


(5.48) 


Step 4. Using 05.4ip . we obtain 


IJ, 


rmain i ^ c 


74+1,72 ^74+1,74! — ^74+1 


+ 1 . n + 1 / h(0) 1 


2n + 1 


rmain ^ r 
^74+1,74 — ^74+1 


2n + 1 \y/^y/n 


+ In+l-r 


Note that 


1 I 2n + l\ 1 

7S V j Ts' 
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Summarising the last two inequalities, we have 




2n + 1 


n 


Steps 1-4 yield Proposition 15.61 (part 2), if uq = no{E) satishes the following. 
Requirement 4. For n > no, the following holds: 

Ch -^h(O) > Ih—^ H— i=h{tf) —h CH2dn 

\/n ^/^ 'n 




2n + 1 


Note that, according to fl5.47|) and Proposition 14.41 Ch — h(0)/\/^ > 0. 


6 Asymptotics for Vun{tn) 

We consider the sequence tk {k G Z) given by Theorem 15.41 and the quantities 

V Vn{tn) = -2cn - c - hi ^ Vyn-k{tn-tk), n>no + l, (6.1) 

\k\<n—\ 


where Vnif) is given by (I5.19p . 

Remark 6.1. Under the assumptions of Remark [5.51 we have Vvnitn) = Vunitn)- 
In this section, we will prove that Vvn{tn) < 0. 

Theorem 6.2. There exists A^j > 0 depending on hi, c, U such that, for all n > uq + 1, 


VUn(tn) + 


3hi 

T” 


< A^n 


^'^nitn) — 


3hi 

~Y' 


6.1 Preliminaries 


Set 

Consider a constant Khi such that 


k 

^n,k • * 

n 


\k\<n—l 


1 1 


h' 



< Khl, 


n > N. 


( 6 . 2 ) 


(6.3) 


Such a constant exists because the left-hand side in fl6.3p is the Riemann sum of a hnite 
integral (note that h'(0) = 0). 
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6.2 Leading order terms 

Substituting tk given by Theorem 15.41 into (16.Ih . we have 


Vvnitn) =-2cn - c-hi ^ Vyn-k{a{n^ - k^) + qn - qk)- (6.4) 

\k\<n—l 

Due to the Taylor expansion, 

Vyn-k{a{n‘^ - k"^) + qn-qk)= - k"^)) + (6.5) 

|fc|<n—1 |fc|<n—1 

where 

:= ^yn-k{a{n^-k‘^)){qn-qk) + ]^ ^ yyn-k{a{n^-k‘^) + in,k){qn-qkf {^-Q) 

|/i:|<n—1 |fc|<n—1 

with < \qn — qk\ < 2En^^‘^. Using f|4.5|] and the functions G{x) = G{a,x) and H{x) = 
H{a ,x) given by fl4.10|) and fl4.11|) . we represent the sum in 06.51) as follows: 

y] Vyn-k(a(n^ - k^)) = n -^i^n,k) + ^ ^H{xn,k) + ^2,^, (6.7) 

\k\<n—l \k\<n—l \k\<n—l 

where 

S2,n := - k,a{n^ - k‘^)). (6.8) 

\k\<n—l 

Now we replace the sums in 06.71) by the integrals. Set (recall that G'(l) = —1/2) 

^ = E (6-8) 

\k\<n—l \k\<n—l 

^h,n ■= y2 -H{Xn,k) - Ih, (6.10) 

\k\<n—l 

where Iq = Ig{o) and Ih = Ini^i) are given by 04.12p . Then 06.7p takes the form 

Vyn-k{a{n^ - k^ = 1^71+^ + ^ + + ^ + S2,„. (6.11) 

\k\<n—l 

Combining 06.4p . 06.5p . and 06.lip and using Lemma Will we obtain 

VUn(tn) = ( —2c — hilG)n + - 2 -~ 4" ^ 2 ,n + nTjg^n + 

= —^- hi + S2,n + n'Eg^n 4—• (6.12) 
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6.3 Remainders and proof of Theorem 16.2 

It remains to estimate S 2 ,n, Sg,n, and in (I6.12p . 

, , f2EBiT3/2 o T 2 

Lemma 6.3. < ^EKhxrT^I^ + -+ ^E^Bl -—^ 

Proof. Using fl6.6D . we write where 


^l,n • 'y ^ V^n—fc(n('n k )){Qn Qk}t 


|A:|<n—1 
1 


•” r) y ^ Vn-kipip ) + fn,k){(ln Qk) ■ 


|A:|<n—1 


Using (14.51) . (14.8p . (16.3p . (I5.17p . and the inequality \qn — Qkl < we have 




E 

\k\<n—l 


E 


-h' 


1 /1 — X. 


n,Ai 


2E5in^/2 


(a(u2 - A;2))3/2 


< 2EKhin-^/^ + 


2EBiT^I2 _i 
a3/2 ^ • 


Further, using (14.9p . (I5.13P and the inequalities \fn,k\ E lo'n — 5'fcl < 2ErP^‘^, we have 


E 


4E^B*n 


^ (a“-(n2 - fc2))2 - (a““) 


2E^B*T2 

n . 


|fc|<n— 




(6.13) 


□ 


Lemma 6.4. |E 2 ,.| < E ..^2 _ 


Proof. Lemma 16.41 is a straightforward consequence of (16.8p and (14.8p . 
Remark 6.5. For n > 4, we have 


□ 


E 

|fc|<n—1 


n-1 


n—1 


rP — /c2 


= n 


-2 


+2^^—^+2u ^(2ln(n—l)+2+n ^) —)■ 0. 


k=l 


k=l 


n — k 


Propositions 14.31 litem 1) and 14.41 imply the following. 

Lemma 6.6. There exists Kg > 0 such that for n > N the following holds: |Eg^„| < Kgn~^^‘^. 
Lemma 6.7. There exists Kh > 0 such that for n > N the following holds: \T^h,n\ < Khn~^^‘^. 
Now Theorem 16.21 follows from (I6.12p and Lemmas 16.3116.71 if the following is satished. 
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Requirement 5. For n > Uq, the following holds: 


2EKhin 


-1/2 , / 2ER1T3/2 

V 


+ 


2E^B*T2 


n ^ + 


E 

\k\<n—l 


Bn 


a(jF — 


+ K,+ 


2 


n 


-1/2 <- 


(6.14) 


7 Estimates of Un{t), Un+i{t),... for t ^ {tn-l, tn) 

7.1 Uniqueness of a switching moment 

As before, we consider the sequence and the functions Vn given by Theorem 15.41 In this 
section, we will prove the following result. 

Theorem 7.1. For all n>nQ + l, we have 

Vn{t)<Q, te[tn-l,tn). (7.1) 

Fix 6q > 0, satisfying the inequality {Kh is given by Lemma [6?7)1 

e^Kh < (7.2) 

The proof of Theorem 17.11 is based on the following proposition. 

Proposition 7.2. Let the assumptions of Theorem 13.21 hold. Then, for all n > no + 1, 

2. Vn{t) < {n - l)"^/2 ^ g [t^_.^,tn-i + 9o{n - 1)^/2], 

2. Unit) > 0 for all t G [tn-i + 0o{n - l)^/^t„]. 

We hrst assume that Proposition 17.21 is true and prove Theorem 17.11 The proof of 
Proposition 17.21 is given in Sections 17.21 and 17.31 below. 

Proof of Theorem U.W By (13.2p (for n = no + 1) and Theorem 16.21 (for n > no + 2), we have 
Vnitn-i) < —3hi/8. Therefore, Proposition 17.21 (part 1) and inequality (17.2p imply that, for 

t e [tn-i,tn-i + 6'o(n - 1)^/^], 

Vn{t) = Vn{tn-i) + ! w(s) ds < + 6'o(n - 1)^/^(n - 1)~^/^ < 0. (7.3) 

Now, for t G [tn-i +6'o(n — 1 )^^^/^]^ Proposition 17.21 fpart 2) imply that Vn{t) can vanish 
no more than once, and thus, by Theorem 15.41 this happens no earlier than at t = □ 

As a corollary of Theorem 17.11 we obtain the following result. 

Theorem 7.3. For all n > no + 1, we have 

Vj{t) <0, t G [tn-i,tn), j > n. 
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Proof. First, we show that 


Vvj{tn) <0, j > n. (7.4) 

To do so, we estimate Avj(t) for t G (tn-i,tn] and j > n. Using (15.19^ . (I2.5p . and the fact 
that 2/n(t) > 0, we have 

Avj{t) = -2c - hi ^ ^yj-k{t - 4) = -2c - hi ^ yj-k{t - tk) < -2c. 

\k\<n—l |A:|<n—1 

In particnlar, Avj(tn) < 0. Together with the relations Vnifn) = 0 (Theorem 15.41) and 
Vn+i{tn) < 0 (Theorem [62]), this yields fl7.4p . 

On the other hand, fl2.6p implies that Vyj-kit — tk) decreases and thns 

Vvjit) =-c{2j+ 1) - hi ^ Vyj-kit-tk) 

\k\<n—l 


increases. Together with (17.41) . this yields Vvj{t) < 0 for all t G (4_i,4] and j > n. Since, 
additionally, Vn(t) < 0 for all t G [tn-i,tn) (Theorem [7TD, the desired resnlt follows. □ 

Remark 7.4. Under the assnmptions of Remark 15.51 Theorem 17.31 implies that 4 is indeed 
the switching moment of the node 

Remark 7.5. Theorem 17.11 implies that the eqnation in (I5.28p has a unique root on the 
interval [—E\/n + 1, U\/ri~+T] 

In the rest part of this section, we will prove Lemma 17.11 

7.2 Proof of Proposition 17.2[ Part 1 

7.2.1 Leading order terms 

We take 6 G [0, Quin — 1)^'^^] and set t = 4-i + 6^ G [4-i, 4-i + 0o{n — 1)^^^]. 

First, we represent Vnitn-i + 0), using (I5.19p and the relation = 4-i, as follows: 

hn(4 —1 T hi 2c hi ^ ^ yn—l—kifn—l T ^ tk) 

\k\<n-2 

-hi 22 ^ yn-l-k{tn-l + 0 - tk) - hi{yi{e) + y 2 n-l{ 0 )) . 

\k\<n—2 


Set 


m \= n — 1. 


Since yi{6) > 0 and y 2 n-i{G) > 0 for all 0 > 0, we obtain 


hn(tm T ^ hi 2c hi 2 ^ ijm—kitm T tk) hiPi^^^^ 

|A:|<m—1 


(7.6) 
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where 

Ss.m := + 6* - 4)- (7.7) 

\k\<7n—l 

Further, to apply the Taylor expansion in fl7.6p . we note that 

tm-tk + 9 = a{m^ - k"^) + (g^ -qk + O)- 

Therefore, using (14.4p . the function H{x) = H{a^x) given by (14.lip , and (l6.1Up . we obtain 


\k\<m—l \k\<m—l 

^ ^ H {^Xffi T S5 T S4 I}j T T S5 ^ T ^4,7; 


\k\<m—l 


where 


^ 4 ,m • ^ ^ ym—ki^^i^^ ) T ^m,/c) (o^m. T , 

\k\<m—l 

^ 5 ,m-= y] ri(m - fc,a(m^ -/c^)), 

|fc|<m—1 

Xm,k = m/k, and —2Em}/‘^ < ^ra,k < 2Em^/‘^ + . 

Combining (17.61) . (17.81) . and Lemmas 14.21 and WER we have 

Vn{tm + 9) <hi {Khfn + ISs^ml + |S4_m| + . 

7.2.2 Remainders 

It remains to estimate Fl4 „j, and ^ in fl7.11|) . 

Lemma 7.6. iSs^^l < Khim~^ + ( + {2E + 9o)Bl , I m“42 




Proof. Using fl7.7p . we write where 

^'z,m ■= 51 Vy„,-k{a{m^ - fc^)), 

\k\<m—l 

^3,m := 51 ^ym-k{a{wf - k^) P ira,k){Qm- qk + 9). 

\k\<m—l 

Using fld.Sp . fld.Sp . fl6.3p . and (15.17p . we have 


|V,„| < 



(7.8) 


(7.9) 

(7.10) 


(7.11) 
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Further, using fl4.9p . fl5.17p . and the inequalities \qm — (Ik + < (2-F + and 

a{m? — + ^rn,k > — k'^), we have 


IS" I < 

I 3,ml — 


Lemma 7.7. ISa™! < 


S 


(2E + eo)m^/^B* , . ^ * T 2 0/2 

^^ < {2E + 9o)B ;-— 


{2E + 6'o)v42T3/2 




□ 


-m 


■ ■ ■ 

Proof. Using fl7.9p . fl4.9p . fl5.17p . and the inequalities \qm — qk + 9\ < {2E + and 

a{m? — k"^) + ^m,k > — /c^), we have 


-'4,m| ^ 


s 

\k\<m— 


(2E + 9o)A2m^^‘^ ^ (2F^ + 6'o)A2^3/2__i 


, (a”-(m2 - fc2))3/2 - (a”-)3/2 


m 


Lemma 7.8. jSs^l < ^ 

Proof. The assertion is a straightforward consequence of fl7.10p . fl4.7p . and fl5.17l) . 


□ 


□ 


Now inequality fl7.1ip together with Lemmas I7.6117.8I yield part 1 in Proposition 17.21 if 
the following requirement is satished 

Requirement 6. For n > uq, the following holds: 


Khi + 


{2E + 6'o)v 42T3/2^ -i- ^o)-B2T2 ^ 1 ^ 3 / 2 ^ ^_i ^ Kh 


(a“-)3/2 


7.3/2 


+ 




+ ^ I ™ ‘ < 


2 

(7.12) 


7.3 Proof of Proposition 17.2[ Part 2 

7.3.1 Preliminaries 

We introduce the function 




+ h"\^ 


1 I + X 


Note that 


Fix 


a V 1 — 

limT(a:) = h"(0) = 

x^l 


X 6 (0,1). 


)) € 0 , 


4v^ 


(7.13) 

(7.14) 

(7.15) 


We will need the following lemma. 
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Lemma 7.9. There exist Xq G (0,1) and £o > 0 such that 


h" 


1 — X 


■\/a{l - x^) + El 
for all X G [xq, (m — l)/m] and ei with 


+ h" 


1 + X 


A/a(r^^x2y+^ 


< —rj 


(7.16) 




(7.17) 


Proof. Choose rji G (?], l/{A^/n)). Equation fl7.14p implies that there is Xq G [0,1) such that 


4'(x) <- 771 , xg[xo, 1 ). 

It is not difficult to check that there exists £0 > 0 such that 


1 — X 


^/a{l — x^) + e 


h" 


1 — X 


\/a(l -2:2) 


< 


Vi-V 


for all X G (—1,1) and £ satisfying |e:| < £0 and £ > —a(l — x‘^)/2. 
Formulas 07.181) and 07.19P imply Lemma 17.91 

We fix xo and £0 from Lemma 17.91 and h such that 

b > 2a. 

We introduce numbers Ri and R 2 satisfying 

1 

(m? - 

\k\<xom 

1 


(m2 — Ipy’R 


-2 


E 


kG[xQmym—l] 


Set 


(a™®'^(m2 — k'^) + bm)^R 


3;S[0,oo) 


(7.18) 


(7.19) 


□ 



(7.20) 

, m > N, 

(7,21) 

m > max ( N, -] . 

(7.22) 

V 1 - 2:0 y 

sup |h""(x)| = h""(0). 

S[0,oo) 

(7.23) 


7.3.2 Leading order terms 

As before, we assume that m = n — 1. We take 9 G [ 6 'om^/ 2 ^ 577 ^] and set t = tm + 9. Then 
t E [tm + 6*om^/2^ and the latter interval contains [tm + 9om}/‘^, R], if the following 

is satisfied. 

Requirement 7. For n > no, the following holds {see O7.20p ): 

a{n — 1)2 — E{n — 1)^/2 _|_ 5(^,7 — 1) > an^ + En^^'^. (7.24) 
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(7.25) 


First, we represent Vn{tn-i + G), using fl5.19p . as follows (cf. fl7.5p i: 

Vnitm “1“ -\- l2,m 

where 

Il,m ■ ^ ^ ym—kitm & ^fc); ^2,m ■ "F i/2m+l(^)) 

\k\<m—l 

^6,m •” ^ ^ ^ijm—k(tm ~\~ d tk)- 

\k\<m—l 

Using fl4.4l) . we represent as follows: 


(7.26) 

(7.27) 


U.m 


E 

\k\<m— 


h" 


m — k 


^ (a(m2 - k^) + qm-qk + 0 )^/^ ^ ^a{m^ - k^) + qm - qk + 0 


where 


“h ^7,my 

(7.28) 

(7.29) 


Sy.m := ^ r2(m - k, a{m^ - fc^) + qm.- qk + 0)- 

\k\<m—\ 

Now we split the sum in 07.281) into two sums in which the summation is taken over 
\k\ < Xom and \k\ E [xom,m — 1], respectively. Let us estimate the first sum, using 07.2ip . 
07.23p . and the inequality a{m‘^ — k‘^) + qm — qk > a™“(m^ — fc^): 


\k\<XQm 


m — k 


(a(m2 - + 61)3/2 


< 




B 


h2 


, Bh2Rl -2 
< -— m . 


(7.30) 


(a™"((m2 - A;2))3/2 - (a““)3/2 


\k\<xom 

To estimate the second sum (which we do if m — 1 > xom, i.e., m > 1/(1 — xq)), we set 

qm-qk + 0 


^m,k • 




Below we assume that the following holds. 

Requirement 8. For n > max (? 7 ,o(F^), ); following holds: 


, a { ( 0 o- 2 E)mi /2 

- mm £o, 7 ^ 2 -< - R -. 

2 m \ mil 


Requirement 9. For n > max (no(i?), the following holds: 


2Em^R + bm 


m^ 


< eo- 
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Then 


2m 


m 


min £o, 7 ^ 2 -< em,k < £o, < 


a(m^ - + Qm- qk 


< a. 


max 

m 


(7.31) 


Using (17.31|1 as well as (I7.16p and (\7.22h . we have (xm,k is given by 


E 


m — k 


\k\^[xQTn^m—l] 


_ (a(m2 - fc2) + + 0)3/2 I ^0(^2 _ ^2) 


E 


kG[xom^m—l] 

h" 


(a(m2 - fc2) + + 0)3/2 


1 — X. 


m,k 


sj^m,k) £m,k ^ 


h" 


1 + X. 


m,k 


< —T] 


E 


1 


(a™'“(m2 — k?) + 


\J^m,k) “I" £m,k ^ 
< —r]R2m~^^^. 


(7.32) 


k£[xom,m—l] 

Thus, using fl7.28p . fl7.30p . and fl7.32p . we obtain the following estimate for R m in fl7.26p : 

h,m < + E 7 ,™. (7.33) 

In what follows, we assume that the following is satished. 

Requirement 10. For n > no, the following holds: 

To < Oon^^"^. 

Now we represent l 2 ,m in (I7.26p . using (14.4p and the equalities h"(0) = —l/(4y'¥) and 
h'"( 0 ) = 0 , as follows: 


^2,m - T^h' ^ + 


:h'‘ 


03/2 V\/0/ 0"/' 

1 ^""(0 un 

H-+ h” 


2m + 1 

~7r 

2m + 1 

~7r 


+ r2(l, 0) + r2(2m + 1, 0) 

+ 03/2(r2(l, 0) + r2(2m + 1, 0)) 


03/2 40F ' 20 

where f G [O,0“^/^]. Below we assume that the following holds. 

Requirement 11. For n > no and m = n — 1, the following holds: 

hdfiA —1/2 I 7// / \ I 27 I 2 —1/2 ^ ^ 

—m ' + sup h [xj H —' < 


20 , 


->2m+l 

- vrT 


0 , 


40r' 


Hence, 

l2,m < 0. 

Due to (I7.25p . (I7.33p . and (I7.34p . Proposition 17.21 (part 2) follows if 


— r]R2m 3/^ + 


Bh2Rl 


m ^ + |S 7 .m| + lEfi.ml E 0 . 


1)3/2 


-‘7,m 


6,m I _ 


(7.34) 

(7.35) 
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7.3.3 Remainders 

Let us prove (17.351) . To do so, we need to estimate S 6 ,m and 

Lemma 7.10. iSgmi ^ 


T^*rp 

±>2-^2 _2 
< ^—r^m . 


7 min'^ 2 
‘'m 


Proof. Using fl7.27p . fl4.9p . fl5.17p . and fl7.3ip . we have 


I < 




BX 


1^2 J -2 -2 

< ^ . 


\k\<m 


( a ““( m 2 - / c 2))2 - ( a “-)2 


Lemma 7.11. IS 


7,m I ^ 




2^5/2 _n 

m . 


(a-“)5/2 

Proof. Using fl7.29p . fl4.7p . fl5.17p . and fl7.3ip . we have 


1^7,771 I < 


E 


iin(^2 _ ^2^ _j_ ^^^l/2^5/2 


A2T5/2 _ 5/2 

< ^— . ,, m ' 


iin^5/2 


Using Lemmas 17.101 and 17.111 we see that fl7.35p holds, if the following is satisfied. 
Requirement 12. For n > no and m = n — 1, the following holds: 


Bh 2 Ri ^ B2T2 


1)3/2 


-,min'\2 


ni ^ < r]R2. 


at 


1)5/2 


(7.36) 


8 Main result: proof of Theorem 13.2 


For n = no + 1, Theorems 15.41 and 17.31 imply that the node Unit) achieves the threshold 0 at 
a time moment tn = an^ + qn, where Qn G [—Ey/n, 0] with the same E as in fl3.ip . Moreover, 
neither of the nodes n 77 ,n 77 +i,... switches on the interval [^ 77 - 1 ,^ 77 ), and thus Unit) switches 
exactly at the moment tn. Furthermore, by Theorem 16.21 the required estimates for Vun{tn) 
hold. Thus, the assertion of Theorem 13.21 holds for n = no + 1. 

In particular, we see that items 1 and 2 in Definition 13. ll hold with no replaced by no + 1 
and with the same E as before. Hence, we can repeat the above argument to obtain the 
assertion of Theorem 13.21 for n = no + 2, and so on, by induction, for any n > no + 1. 


A Equivalence of three equations: proof of Proposi¬ 
tion 14.21 

We prove that 

1. equation fl4.15p has a unique root. 
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2. equations fl4.13p and fld.lSp have the same roots, 

3. equations fl4.14p and fl4.15p have the same roots. 


We will prove in detail items 1 and 2. 

Making the change of variables y = —;=W-in fl4.11|] and fl4.10l] . we have 

mV 1 + X 


-fir (a) — 


2h{y) 

lo 1 + 


dy, Ida) = 


'■ + 0O 


4:ayg{y) 

lo (l + ay^)^ 


dy, Ipia) = 


f+oo q^2„,2 


8a y f{y) 
lo (l + ai/2)^ 


dy. 


Now we see that Inia) decreases from 1 to 0 as a increases from 0 to +cx). Hence, for any 
0 < c < hi/2, equation fl4.15p has a unique root a > 0. Item 1 is proved. 

Let us prove item 2. Integrating by parts and using equations fl4.1|) and fl4.2|) . we obtain 


'‘+00 


Ipia) = 


2,,2 


8ad 


(1 + ay"^)^ 1 + ay"^ 


f{y)dy + 


<■+00 


'0 


1 + ay 


:f{y)dy 


ay{-l + ay"^) 

(1 + 


yiy) 


H-cxd 


'*+00 


ay{-l + ai/2) 


ay(-i + ay^) 

^^^9{y)dy + 

/• + CXD 


^+OD 


(1 + ay^ 


-J{y)dy 


+ y- 


a 


(1 + ay"^)"^ 1 + ay"^ 

+ CX3 /*+00 


9{y)dy + 


a 


'0 


1 + ay 
1 

'2 


-.9{y) 


+ 


1 


My)dy + 


1 + ay'^ 
2a 


1 + ay‘ 

{f{y) -y9{y))dy 


0 


1 + ay'^ JO 

n = ‘^^^Hia) -1- 

lo 1 + ay^ 2 2 


1 + ay^ 


h{y)dy 


It is easy to conclude from the last identity that equations (14.131) and fl4.15p have the 
same roots. Item 2 is proved. 

Similarly to item 2, integrating by parts and using relations (14.ip and (14.2p . we obtain 


0+00 


Ida) = 


Aay 


lo (1 + 02/2)2 


9{y)dy = 


1 + ay- 


^9{y) 


+ 00 


^+CO 


+ 


+ 1 + ay 


My) = -1 + Ida). 


From this identity, it is easy to conclude item 3. 


B Requirements on no(^) 

In this appendix, we collect the constants that we use throughout the paper to determine 
no{E) as well as all the 12 requirements on the number no = no{E) entering Dehnition 13.11 
of admissible E. 
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B.l Constants not depending on a or E 

1. To > 0 is an arbitrarily fixed real number (see (I4.6p ). 

2. Set (see (ED) := sup„>o,i>^o Vt yn{t) - Vi f 

3. Set (see (112D) -4i := tVi \^yn{t) - ^ ^ 

4. Set (see (ED) := sup^^o,t>To ^Vt) - Vi f ^ 

5. Set (see (ED) ^ := sup„>o,t>^o jjnit) - 

6 . Set (see (gSD) Bq := sup^>o,t>ro t V?/„(t) - g ^ ^ 

7. Set (see (ED) := sup^>Qt>^jVi ^VnE) - \h' 

8. Set (see (ED) := sup„>o,t>^o t\/t|y„(t)|. 

9. Set (see (ED) B^ := sup„>o,t>^o t^|Vi/n(t)|. 

10. We use the notation (see (15.181) ) R := vr. 


11. Bh 2 is given by (see (17.231) ) Bh 2 '■= sup \h"{x)\ = -h"(0). 

a;S[0,oo) 


12. BhA is given by (see ([TED) Bm ■= sup |h""(x)| = h""(0). 

a;S[0,oo) 

B.2 Constants depending on a but not depending on E 

1. a is a unique root of equation (I4.13p (or equivalently (14.141) . (I4.15p ). 

2. Consider the function ha{x) := h{^ + . Set (see (ED) Da := inCe(o,i] ha{x). 

3. Set (see ([ED) P ■= sup^6(o,i] 

4. iV is a fixed natural number satisfying (see fl5.8p . fl5.9p ) 

A'(m“ -1) -(2^ -1) -(2v^ -1+ n (^ s ^)^ - 1)) (2^ -1) ^ > 0. 

5. Set (see ( 1 ^ ) Dpi ■.= 2V - 1 , Dp2 ■= A^((l + - l). 

6. Constants needed to define K and K' from Lemma 15.31 


(a) Set (see ([E2D) Li{F,N) := ^^^F{x)dx -Y.\k\<n-i^B{^ 

the supremum exists due to Proposition 14 31 Set K[ := hiLi{F, N). 


, where 
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(b) Set (see fl5.33p ) li := hiA^-^R. 

(c) Set (see (|5.34p ) K' := K[ + 

(d) Set Ki ■= sup„> 7 v'n-^/^|C'^+^ — Cn^, where 

Va(/2 -P)/|- 

\ ifci<^-i V 


I - k 


^a{P - k^) 


, / = n, n + 1. 


(e) Consider the function Fix) := , ^ f ( , ^ ^ ) and constants L L*., I? > 0 


Note that the supremum exists due to fl5.35l) . 


such that (see fl5.37p i 
f\/E An- 


j F{x)dx+ Y. 

-n<k<n-l ^ ^ 


n> N. 


Set (see (I5.38P ) K 2 := h (Lj - L*- + ■ 

(f) Set (see (| 5.39p ) K 3 := 2hiii^T3/2. 

(g) Set (see ([ElQ])) K-=K^ + K 2 + K^A. 

Remark B.l. In principle, due to Proposition 15.31 we could dehne K' := sup„>^ 
and K := sup^>^n^/^|Cn+i — Cn\- However, calculation of the values Cn is compu¬ 
tationally consuming as it involves Bessel functions. We used the strategy described 
above, since it is based on error estimates of Riemann sums only. 


7. Set (see flR8l) ) x := Dp 2 - Dpi - 2 {Dpi + Dp 2 )DpiA (> 0). 


8. Set (see (EEl) Bo - (> 0). 

9. For a > 0, set (see l|aja) So := sup„>^(VnS|i|<„-i ““’y 

values of Si, S 2 , and S 3 to determine nQ{E). 

10. For a > 1, set (see (I5.17P ) T„ := sup„>^(Xl|fc|<n-i (na-I^)") • values of 

T 3 / 2 , T 2 , and T 5/2 to determine no{E). 

11. Let constants L^i^H, N) > -^h(O) and iV) > 0 be such that (see Proposition 14. 4p 

1 1 n-l 

li{h,n)^-12{h,n)-<Ih- 

n n 

k=—n 

where H{x) is given by fl4.1ip and In = Inici) is given by fl4.12p . Set (see fl5.47p l 
Ch-.= LI{H,N),Ih-.= 12{H,N). 


n 


n > N, 
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12. Consider the function H{x) := y^H{x). Set (see fl5.48p i 


Ch2 :=sup„>jvn 


13. Set (see 06.31) i Khi := sup„>jv’^ 


-1 


E 

\k\<n— 


^ a(l-(fc/n)2)^ (v^\/l 


—kin 

1+k/n 


14. Set (see Lemma ICHl and (16.Qh ) Kg := sup„>^n^/^ 
G{x) is given by 04.1Up . and Iq = Ig{o) is given by ([ 


|fc|<n—1 

'M- 


, where 


15. Set (see Lemma ICTI and 06.10p i Kh := sup„>jv^^'^^ 
is given by (14.lip , and Ih = Ih{cl) is given by (I4.12f 

16. 6*0 satishes (see 07.2p ) 9oKh < j- 

17. T] satishes (see (17.151) ) rj G ^0, • 


E 

|fc|<n—1 


, where H{x) 


18. xo G [0,1) and £o > 0 satisfy (see 07.16p . fl7.17p ) 

1 — X \ 1 + X 


h" 


for all X G [xq, (n — l)/n] and ei with 


+ h" 


A/a(l - x^) + El 


< —rj 


min ( ^0, ^ (^2 - ^ < £o, n > max (^N, ^ 


19. b satishes (see 07.201) ) b > 2a. 

20. Set (see (|7.2ip) := ^ („2_fc2)3/2 

\k\<xon 

B.3 Constants depending on E 

1. Set (see := a - 

2 Enl/2 


2. Set (see (EH) ■= a + 

3. Set (see 05.15P ) 5n := + 


_Ml! 2A*£ai/^ 1 

aDa n ' Da(a™in)3/2 y/n‘ 


4. Set (see (| 7.22p ) R 2 := inf„>max(v,i/(i-xo)) E 


k(£[xon,n—l] 


_ n3/2 

(^max ^^2 _/j2^_|_J„p/2 ' 
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B.4 Requirements on no(^) 

We assume that the following requirements hold for n > no{E): 
1. n > N. 


9 „min \ TO 
“n — 2n-l‘ 


3. dn + 5n+l < 2Zi)pi. 


4. + ^ + + ;4(l + fJ (M0)|±ii4„-n + (1 + 

5. 2EKarr’i + ( + E a(„f"t2) + (^s + £ §■ 

^ V n ; / |fc|<n-l 

( K I (2^+®o)"^2h3/2 \ -1/2 I f B 1 T 3/2 I (2£;+6»o)B|T2 I Ai73/2\ 1 ^ 

0- + (a--)3/2 J n + + (amin )2 +^;^jn S—■ 

7. a{n — 1)^ — E(n — 1)^/^ + b{n — 1) > an^ + Ev}!'^. 

8. ^ <min(£o,|^(2-i)) for n > max (no(^), . 

9. 2En~‘^/‘^ + bn~^ < Sq for n > max (^no{E), • 

10. To < 6*on^/^. 


11. 1/2 + sup^>2:^ h" (x) + 1/2 < 1 . 


1 9 / ^fe2^1 

• I (a“i")3/2 


I ^2^2 \ -1/2 I ^2T5/2 1 ^ „P 

+ Ja^ + (a--)5/2^ S V^2- 
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